Local Atomic Strain in ZnSci ^Te^ from High Real Space Resolution Neutron Pair 

Distribution Function Measurements 
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High real-space resolution atomic pair distribution functions (PDFs) have been obtained from 
ZnSei_a;Tea; using neutron powder diffraction. Distinct Zn-Se and Zn-Te nearest neighbor (nn) 
bonds, differing in length by Ar = 0.14 A, are resolved in the measured PDF allowing the evolu- 
tion with composition of the individual bond-lengths to be studied. The local bond-lengths change 
much more slowly with doping than the average bond-length obtained crystallographically. The nn 
bond-length distributions are constant with doping but higher-neighbor pair distributions broaden 
significantly indicating that most of the strain from the alloying is accommodated by bond-bending 
forces in the alloy. The PDFs of alloys across the whole doping range are well fit using a model 
based on the Kirkwood potential. The resulting PDFs give excellent agreement with the measured 
PDFs over the entire alloy range with no adjustable parameters. 
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I. INTRODUCTION 

Ternary alloys, such as ZnSei_a;Te2;, are technologi- 
cally useful because the band-gap can be tuned between 
the end-member values as the composition, x, is varied. 
This makes the proper characterizationjpf these materi- 
als the subject of much investigation.Ercl ZnSei_a;Te^ is 
an example of a II- VI semiconductor pscudobinary alloiz: 
that can be made over the entire range of compositions.u 
II- VI alloys are becoming increasingly important because 
they are often used as the basis for magnetic semiconduc- 
tors with the additional alk)ying of small amounts of Mn 
on the metal sublatticcQlZl The recent suggestion that 
high-speed logical circuits can be made out of devices us- 
ing spin diffusion instead of electron diffusion (so-called 
"spintronics" ) is, adding extra impetus to research into 
these materialsJa Clearly, it is important to be able to 
characterize the atomic and electronic structure of these 
alloys in detail. 

The study of alloys is complicated by the fact that 
considerable local atomic strains are present due to the 
disordering effect of the alloying. This means that local 
bond-lengths can differ from those inferred from the sf^-, 
age (crystallographic) structure by as much as 0.1 ABt£l 
This clearly has a significant effect on calculations of elec- 
tronic and transport properties 13 To fully characterize the 
structure of these alloys it is necessary to augment crys- 
tallography with local structural measurements. In the 



past the extended x-ray absorption fiHc^tmcture (XAFS) 
technique has been extensively usedJjOEj More recently 
the atomic pair distribution function (PDF) analysis of 
powder diffraction data has also been applied to get ad- 
ditionaLjkjcaJ structural information from In^jGai-xAs 
alloysJlJlljO In that case high energy x-rays combined 
with good resolution and a wide range of momentum 
transfer allow the In-As and Ga-As nearest neighbor 
peaks to be resolved. In this paper we describe PDF mea- 
surements of the II- VI alloy ZnSei-a^Te^: from neutron 
powder diffraction measurements using the new General 
Materials Diffractometer (GEM) at ISIS. In these mea- 
surements the distinct Zn-Se and Zn-Te bonds, which dif- 
fer in length by just Ar = 0.14 A, could be distinguished 
demonstrating the quality of the data from GEM. 

Both ZnTe and ZnSe have the zinc-blende structure 
(F43m) where the Zn atoms and Te, Se atoms occupy the 
two interpenetrating face-centered-cubic (fee) lattices. In 
the alloys the lattice parameter of ZnSei_a;Tea; inter- 
polates linearly betweeia^he end member values consis- 
tent ■vjdi.k.Vegard's law f^l. However, both XAFS experi- 
mentstirEj and theoryt3~0 show that the atomic nearest 
neighbor (nn) distances deviate strongly from Vegard's 
law. Rather, they stay closer to their natural lengths 
found in the end-member compounds: L",, ^^=2.643(2)A 
and L°„_s^=2.452(2)A. 

A limitation of the XAFS method for studying the lo- 
cal structure of alloys is that it only gives information 



about the first and second neighbor bond-lengths and in- 
formation about the bond-length distributions with less 
accuracy. In this study the PDF analysis of neutron pow- 
der diffraction data is used. The PDF is the instanta- 
neous atomic number density-density correlation func- 
tion which describes the atomic arrangement of the ma- 
terials. It is the Sine Fourier transform of the experimen- 
tally observed total scattering structure function, S{Q), 
obtained from a powder diffraction experiment. Since 
the total scattering structure function includes both the 
Bragg and diffuse scattering, the PDF contains both local 
and average atomic structure yielding accurate informa- 
tion on short and intermediate length-scales. Previous 
high resolution PDF studies on In^^Gai—tjAs were carried 
out using high energy x-ray diffraction.E3o This yielded 
data over a wide Q-range {Q is the magnitude of the 
scattering vector) which resulted in the very high real- 
space resolution required to separate the nearest neighbor 
peaks from In-As and Ga-As. The high Q-range cover- 
age and Q-space resolution of the new General Materials 
(GEM) Diffractometer at the ISIS neutron source allowed 
us, for the first time, to obtain similarly high real space 
resolution PDFs of ZuSei-xTSx using neutrons and to re- 
solve the Zn-Se and Zn-Te bonds that differ in length by 
only 0.14 A. Furthermore, the data collection time was 
only sixty minutes compared to the 12 hours for the x-ray 
data with similar quality. The nn distances and average 
peak widths are fit using model independent techniques. 
The PDFs of the full alloy series have been calculated 
using a model based on the Kirkwood potential giving 
excellent agreement over a wide range of r with no ad- 
justable parameters. 



II. EXPERIMENTAL 
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FIG. 1. (331) and weak (420) peaks of ZnSei_:,Te^ mea- 
sured at 300K using Cu rotating anode. 



peaks are reproduced in Fig. n^. The double-peaked shape 
comes from the Kd and Kq,2 components in the beam. 
The line-widths are narrow and smoothly interpolate in 
position between the positions of the end-members veri- 
fying the homogeneity of the samples. 



A. Synthesis and Characterization 

Finely powdered samples of ~10g of ZnSei_a;Te2: were 
made with x=i, |, |, |, |. The starting reagents (zinc 
selenide, metal basis, 99.995%; zinc telluride, metal ba- 
sis, 99.999%) were finely ground, mixed in the correct 
stoichiometry, and sealed in quartz tubes under vacuum. 
The samples were then heated at 900°C for 12-16 hours.E2l 
This procedure (grinding, vacuum sealing, and heating) 
was repeated four times to obtain high quality homo- 
geneous products. The colors of the solid solutions vary 
gradually from dark red (ZnTe) to yellow (ZnSe) as the x- 
value decreases reflecting the band-gap of the alloy sam- 
ples smoothly changing in the optical frequency range. 
The homogeneity of the samples was checked using x-ray 
diffraction by monitoring the width and line-shape of the 
(400), (331), (420), and (422) Bragg peaks measured on 
a rotating anode Cu K„ source. Finely powdered sam- 
ples were sieved through a 200-mesh sieve then packed 
into flat plates and measured in symmetric reflection ge- 
ometry. The (331) and (weak on high angle side) (420) 



B. Neutron Measurements and Data Processing 

Time of flight neutron powder diffraction data were 
measured on the GEM diffractometer at the ISIS spalla- 
tion neutron source at Rutherford Appleton Laboratory 
in Oxfordshire, UK. The finely powdered ZuSei-^^Te^; 
samples were sealed inside extruded cylindrical vanadium 
containers. These were mounted on the cold-stage of a 
helium cryostat immersed in cold He gas in contact with a 
liquid He reservoir. The temperature of the samples was 
maintained at lOK using a heater attached to the cold- 
stage adjacent to the sample. The empty cryostat, an 
empty container mounted on the cryostat and the empty 
instrument were all measured, allowing us to assess and 
subtract instrumental backgrounds. The scattering from 
a vanadium rod was also measured to allow the data to 
be normalized for the incident spectrum and detector ef- 
ficiencies. Standardpdata corrections were carried out aa 
described elsewhereEllEj using the program PDFgetN.o 
After being corrected the data arc normalized by the total 
scattering cross-section of the sample to yield the total 
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FIG. 2. Q[S{Q) - 1] for ZnSei-;^Te^ measured at lOK. 
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FIG. 3. G{r) — 4TTr{p{r) — po) for ZuSei-xTe^ measured 
at lOK. 



scattering structure function, S{Q). This is then con- 
verted to the PDF, G{r), by a Sine Fourier transform 
according to 



G{r) 



Q[S{Q)-l]sm{Qr)dQ 



47rr(p(r) - po), 



(1) 



where p{r) is the microscopic pair density of the sample 
and Po is the average number density. 

The GEM instrument yields useful diffraction informa- 
tion up to a maximum Q of greater than 90 A~^. Unfor- 
tunately, due to a neutron resonance in Te we were forced 
to terminate the Fourier transform at a maximum Q of 
40 A"'^ in this experiment. This resulted in nn peaks in 
these alloys which are resolution limited rather than sam- 
ple limited. This was verified by Fourier transforming the 
ZnSe end- member at higher values of Qmax ■ The nn peak 
kept getting sharper up to Qmax — 60 A~^. Nonethe- 
less, the distinct short and long bond distances are still 
evident in the alloy PDFs. At the time of this measure- 
ment the backscattering detector banks on GEM were 
not operational. With the backscattering banks yielding 
better statistics in high Q and adding detector cover- 
age, one might expect to get similar quality PDFs in a 
fraction of the time. The reduced structure functions, 
Q[S{Q) — I], obtained from the ZuSei-^Tca: samples are 
shown in Fig. and the resulting PDFs are shown in 
Fig. I 



C. Method of Modeling 

The PDF can be calculated from a structural model by 
taking advantagp-pf the definition of radial distribution 
function (RDF),Eil T{r), 



T{r) ~ Anr p(r) 






,), 



(2) 



and substituting the calculated p{r) into Eq. |[ Here bi 
is the scattering length of the zth atom, (6) is the scat- 
tering length averaged over the sample composition and 
isotopes, rij = jr^ — rj\ is the distance separating the zth 
and jth atoms and the sums are taken over all the atoms 
in the entire sample. To approximate the atomic thermal 
motion we convolute the delta- functions with Gaussians. 
Before being compared to the data the calculated G{r) 
is convoluted with a termination function, si'ii{Qmaxf)/fi 
to account iot-the effects of the finite range of the mea- 
sured data.cj'Ej 

We have used three approaches to obtain structural 
information from the PDF. First, we carry out a model 
independent analysis by fitting Gaussian functions to 
peaks in the RDF. Next we calculate the PDF expected 
from the average crystal structure and refine, using a 
least squares approach, atomic displacement (thermal) 
parameters to obtain empirically the PDF peak widths 
in the alloys. This-, is done using the PDF refinement 
program PDFFIT.E£l Finally, we calculate the PDF from 
a Kirkwood potential based model where the atom po- 
sitions and thermal broadenings are fully determined by 
the atomic potential parameters. 



The ZnSe and ZnTe nn distances in the alloys were 
found by fitting two Gaussians, convoluted with termi- 
nation functions to account for the termination effects, 
to the nearest neighbor peak in the PDF. Peak positions 
and widths were varied. The relative peak intensities 
were constrained to those expected from the alloy com- 
position. The widths refined to the same values as the 
end-members within the errors for all the alloys. We thus 
repeated the fits constraining the peak widths to have the 
values refined from the end-member PDFs. This more 
highly constrained fitting procedure resulted in less scat- 
ter in the refined peak positions. i— . 

To find the far neighbor peak widths PDFFITnJ was 
used. The zinc-blende crystal structure was used with 
all the atoms constrained to lie on their average posi- 
tions. Lattice parameters, scale factor, isotropic thermal 
facto rs ,, a nd r-dependent PDF peak broadening param- 
etersESE3 were allowed to vary but the atoms were not 
allowed to move off their sites. In this way, all of the 
atomic disorder, static and dynamic, is included in the 
refined thermal factors that are giving an empirical mea- 
sure of the PDF peak widths at higher-r. This approach 
is better than fitting unconstrained Gaussians because 
of the problem of PDF peak overlap at higher distances. 
Clearly this approach does not result in a good model 
for the alloy structure but will yield a good fit to the 
intermediate PDF in the alloys. However, it should give 
reliable empirical estimates of the width of high-r PDF 
peaks, even when they are strongly overlapped, and al- 
lows us to separate the disorder on the cation and anion 
sublattices. The PDFs were fit over the range of 3 to 
15 A~^. This range was selected so the global properties 
of the alloys could be fit without influence from the nn 
behavior. 

Potential based modeling to yield realistic alloy struc- 
tures has been carried out using a model based on the 
Kirkwood potentialjtl-I'lii.'Mjrocedure has been described 
in detail elsewhere. ll30c3'E3 The model consists of 512 
atoms arranged in the zinc-blende structure with peri- 
odic boundary conditions where the interatomic force is 
described by the Kirkwood potential. The system is then 
relaxed by moving atoms to minimize the energy. 

The Kirkwood potential can be written as. 






^L'.,? + ^LlJ2 (cos ^ 
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V 
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where L^ is the nn bond length of an undistorted ref- 
erence crystal structure, Lij is the length of the bond 
between the atoms i and 7, and L?, is the natural bond- 
length. In this definition the bond-stretching, a, and 
bond-bending, /?, force constants have the same units 
and 9ijk — am:os(— -1), for an ideal tetrahedron. Lit- 
erature valuesLj were used for the bond-stretching and 
bond-bending parameters obtained from elastic constant 
measurements. We also tried optimizing a and f3 by fit- 
ting to the ZnSe and ZnTe end-member PDFs; however, 
the PDFs calculated using both sets of parameters gave 
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for 


ZnSei- 


^Te^ and 




a (N/m) 




(3 (N/m) 




P/a 


a** 


ZnSe 
ZnTe 


33.7 
31.1 




4.6 

4.7 




0.14 
0.15 


0.78 
0.76 


InAs 
GaAs 


35.1 
44.3 




5.8 
9.2 




0.16 
0.21 


0.74 
0.70 



comparable agreement when compared to the alloy data 
so we simply report the results obtained with the litera- 
ture values of a and /3. The values of a and /? for ZnTe 
and ZnSe used are shown in Table |[ 

The PDFs for the alloys are calculated with no ad- 
justable parameters using the same potential parameters 
used for the end-members. The additional bond-bending 
parameters present in the alloy due to Te-Zn-Se type con- 
figurations are determined as a geometric mean of the (3 
parameters for the end-members. E3 The thermal broad- 
ening of the PDF is calculated by determining the dy- 
namical matrix from the potential and projecting out the 
atomic displacement amplitudes for each phonon.ES 



III. RESULTS AND DISCUSSION 
A. Model Independent Results 

Upon inspection of the PDFs presented in Fig. ^ 
one will immediately notice the splitting of the first 
peak. This comes from the fact that the nn distances 
of ZnTe and ZnSe stay close to the end-member values 
of 2.643(2)A and 2.452(2)A respectively. The positions 
of each component of the doublet were determined by 
fitting Gaussians as described above. The values for the 
nn bond-lengths are shown as filled circles with 2a er- 
ror bars in Fig. H. Also plotted in the same Figure as 
open circles are the nn bond-lengths determinedJrom an 
earlier XAFS study by Boyce and MikkelsonJl3 There 
is clearly excellent agreement between the two results. 
Superimposed on the data are lines which are the pre- 
dictions of the Kirkwood model for the nearest neighbor 
bondplengths using the potential parameters given in Ta- 
ble |.tZl Again, there is excellent agreement with the data 
with no adjustable parameters. 

In contrast to the local structure, the long-range struc- 
ture is well jdeaeribed by the virtual crystal approxima- 
tion (VCA)EI'E3 The VGA assumes that the structural 
properties of a crystal alloy all are a linear interpolation 
of the end-member values. If this were true in semicon- 
ductor alloys then not only would one be able to find the 
lattice parameter of the alloy from Vegard's law but the 
nn distance would be, for zinc-blende crystals Ai_2.B2:G, 

Lac = Lbc — ^o-- This is shown in Fig. as the 
dashed line. 

It is clear that the local bond-lengths remain closer 
to those in the end-members than to the prediction of 




FIG. 4. nn positions from the PDF|4|.), XAFS dataB (o), 
and Kirkwood model z-plot (solid line)LJ as a function of com- 
position, X, for ZnSei-ajTe^,. The dashed line is the average 
nn distance. Note that not all XAFS points had reported 
error bars so they were all set to the same value. 



Vegard's laW|-Jn fact the bond-lengths stay close to the 
Pauling limitcJ in which they would remain completely 
unchanged in length across the alloy series. The devia- 
tion from the Pauling limit is attributed to the disorder 
in the force constants as we describe below. 

The difference in the local structure of the alloys leads 
to a large amount of atomic strain resulting in much 
broader PDF peaks in the high-?- region of the PDF. 
In Fig. H the strain is quantified. The dotted line at 
(Tp =0.0056A^ represents the mean square width of the 
PDF peak, Up, attributed to the thermal motion of 
the atoms while the parabolas indicate additional peak 
broadening due entirely to static strain in the system. 
The largest peak broadening is seen in the mean-square 
width of the Zn-Zn peaks (cation-cation) which is as 
much as 5x as large as the mean-square width due to 
thermal broadening and zero-point motion at 10 K. It is 
also evident that the disorder is larger on the unalloyed 
(Zn) sublattice than the mixed (ScpTei sublattice similar 
to the observation in In2:Gai_2;As.E3'Eil 



B. Kirkwood Model 

The Kirkwood potentialnZl is widely used to describe 
semiconductor alloys. Petkov et ali3 showed that the 
Kirkwood model is good at describing Ina;Gai_a:As, a 
III-V semiconductor. ZnSei-xTex is a II- VI semiconduc- 
tor so it is of great interest to know whether or not the 
Kirkwood model is equally successful for this more polar 
semiconductor alloy. The values of a and /? used in this 
study are shown in Table m with with appropriate values 
for Ina;Gai_a:As for comparison. 
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FIG. 5. Square of the PDF peak widths for far neighbors 
as a function of composition, x. The extracted values (points) 
are plotted with parabolas (lines) to guide the eye. The dotted 
line is at 0.0056A^. 



Other quantities of note are tha-jratio, /3/a, and the 
topological rigidity parameter, a**, t^ which is a function 
of P/a: 



a = 



1 + 1.25(/3/a) 



l + 3.6(/3/a)-f 1.17(/?/a)2 



(4) 



The topological rigidity parameter, a**, can vary from 
to 1 and quantifies the effect of the lattice. The Paul- 
ing limit results when a** = l, the floppy lattice limit. If 
a**=0 the lattice is rigid and Vegard's law will hold true 
locally as well as globally. As can be seen in Table |, the 
values found for a** are close to 0.75 which appears to 
be fairly universal for all semiconductors.Ej 

In this study, the z-plot shown in Fig. H and PDFs 
from each alloy composition were all calculated using the 
Kirkwood parameters appropriate for the end-members 
with no adjustable parameters. Fig. || shows PDFs ob- 
tained from the Kirkwood model plotted with measured 
PDFs for characteristic compositions. The model is very 
successful in matching both the short and longer-range 
behavior of the PDFs for all alloy compositions. The 
measured and calculated PDF peaks of the nearest neigh- 
bor bonds are shown on an expanded scale in Fig. H. It 
is clear that the model based on the Kirkwood potential 
does a very satisfactory job of explaining both the PDF 
peak positions and widths and appears to produce a very 
satisfactory model for the structure of these II- VI alloys. 



C. Comparison wth IriajGai j^As 

The results found for ZnSei_a;Tea; are not entirely 
unexpected. In a previous-, study of Ina;Gai_a;As sim- 
ilar results were obtained.E3 With high real space res- 
olution measurements now possible, direct observation 
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FIG. 6. Comparison of the Kirkwood model (lines) and 
data (o) PDFs for (from top to bottom) ZnTe, ZnSe3/gTe3/5, 
ZnSe4/6Te2/6, and ZnSe. 




FIG. 8. Comparison of the theoretically calculated 
ZnSei-i^Tei: (solid line) and Ina;Gai_2,As (dashed line) 
z-plots. 
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tetrahedral covalent bonds 



literature when atoms are in 



A-B 


rA(A) 


re (A) 


Ta+b (A) 


GaAs 


1.26 


1.18 


2.44 


InAs 


1.44 


1.18 


2.62 


ZnSe 


1.31 


1.14 


2.45 


ZnTe 


1.31 


1.32 


2.63 




FIG. 7. Comparison of the Kirkwood model (lines) and 
data (points) nn distances for ZnSei-aiTej; where a; is (x), 



(O), I (e), I (*), I (A), and 1 (- 



of the nn distances in addition to the static strain in 
the system is observed. The basis for the comparison 
between ZnSei_j;Tej; and In^Gai-ajAs is that the two 
systems are both semiconductor alloys with zinc-blende 
structures. However, they vary in a couple of important 
aspects. The salient difference is number of valence elec- 
trons. ZnSe i_ a; Tea: are more polar alloys and so might 
be expected to have bonding with more ionic character 
than Ina;Gai_2:As which should give rise to smaller (3 val- 
ues since ionic bonding is less directional than covalent 
bonding. Indeed the bond-bending magnitudes are less 
in ZnSei_a;Tea; than In2;Gai_a;As (Tab. |). However, the 
nearest neighbor bonds are stiffer in Ina;Gai_a;As. This 
is presumably also due to the lower polarity of this ma- 
terial resulting in greater orbital overlap and covalency. 
The result is that the (3 /a ratio and a** are similar for the 
two systems. Since it is this ratio, rather than the values 
of a and (3 themselves, that has the greatest impact on 
the structure of the alloy, we find very strong similarities 
in the local structures of Ina;Gai_a;As and ZnSei_a;Tej;. 
For example, the z-plots from both systems are plotted 
on the same scale in Fig. pi. 

The similarity in the z-plots is even more striking be- 
cause of the similar bond-lengths of the end-members in 
these two alloy systems. The ionic radii for the atoms in 
these two alloy series are reproduced in Tabic p. Despite 
the ionic radii themselves being somewhat different, it is 



clear that the sums of the ionic radii of the end-members 
yield values that are within 0.01 A of each other. 

As well as the z-plots of Ina;Gai_xAs and ZuSei-^^Te^^ 
matching rather well, the observed magnitude of the 
mean-square width of the high-r PDFs peak are rather 
similar in the two alloy series. This can be seen by com- 
paring Fig. ^ with Fig. 4 in Rcf. |l^. For example, the 
static strain contribution to the PDF peak widths of the 
unalloyed site has a maximum at 0.027A^ and 0.023A^ 
for X = 0.5 in ZuSei-^^Te^; and In^^Gai-a^As, respectively. 

It thus appears that the structure of the alloys is prin- 
cipally determined by the differences in bond-length of 
the end-members and by the P/a ratio rather than by 
the absolute values of a and (3 or the absolutej-vakies of 
the ionic radii themselves. It has been showrOO that 
the P/a ratios of tetrahedral semiconductors are some- 
what universal resulting in a** values close to 0.75 for a 
wide range of semiconductors. It also does not appear to 
matter whether the cation or anion sublattice is alloyed. 
The unalloyed sublattice accommodates the majority of 
the atomic scale strain. 



IV. CONCLUSION 

From high real space resolution PDFs of ZuSei-a^Te^: 
we conclude the following. In agreement with earlier 
XAFS results and the Kirkwood model the Zn-Se and Zn- 
Te bond-lengths do not take a compositionally averaged 
length but remain close to their natural lengths. Direct 
measurement of this was allowed by the new GEM instru- 
ment at ISIS. The bond- length mismatch creates consid- 
erable local disorder which manifests itself as broaden- 
ing in the PDF peak widths and can be separated into 
thermal motion and static strain. ZuSei-^TCa; was com- 
pared with In^Gai_a;As. Despite having different polar- 
ity the atomic strains in both systems are very similar 
and both are well modeled by the Kirkwood potential 
based model. This suggests that the atomic strains in 
tetrahedral semiconductor alloys are quite universal de- 
pending principally on the bond-length mismatch of the 
end-members and the ratio of the bond-stretch to bond- 
bending forces. 
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